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fSJ ' Abstract. We analyse a high resolution simulation of the formation of a 

cluster's dark-matter halo in a ACDM cosmology (Springel et al. 2001). 
The resolution achieved allows us to map the phase-space structure in 
detail, and characterize its evolution and degree of lumpiness. Scaling 
l/^ ' down the cluster halo to a Milky- Way size halo, we probe the substructure 

CN i expected in the solar neighbourhood. Here we specifically address the 

relevance of such substructure for direct detection experiments aimed at 
^ , determining the nature of dark-matter. 

, 1. Introduction 

in 

I Over the last twenty years a theory has emerged for the formation of structure 

in the Universe (Peebles 1974; White &: Rees 1978). The hierachical paradigm 
has allowed astronomers to make very definite predictions for the properties of 
galaxies today and about their evolution from high redshift. Direct comparisons 
to observations have shown that this model is quite successful in reproducing 
both the local and the distant Universe. However, it relies on a basic assumption: 
2 . that most of the matter in the Universe is dark, in the form of some yet to be 

^ ' identified weakly interacting elementary particles. 

. Thus the crucial test of this paradigm consists in the determination of the 

nature of dark-matter through direct detection experiments. Among the most 
promising candidates from the particle physics perspective are axions and neu- 
^ I tralinos. Axions can be detected through their conversion to electric photons in 

. the presence of a strong magnetic field. The most important direct detection 

process of neutralinos is elastic scattering on nuclei, and the idea is to determine 
the count rate over recoil energy above a given (detector) background level. The 
experimental situation has been improving rapidly over the past years, as the 
relevant regions of parameter space for the different dark-matter candidates are 
starting to be probed (Bergstrom 2000). In all experiments, the count rate is 
strongly dependent on the velocity distribution of the incident particles. In most 
cases, an isotropic Maxwellian distribution has been assumed (e.g. Preese, Prie- 
man & Gould 1988), although there are other examples in the recent literature, 
discussing multivariate Gaussian distributions (e.g. Green 2000). Attempts at 
understanding the effect of substructure on the velocity distribution of dark- 
matter particles have also been made (e.g. Hogan 2001). However, most of this 
work did not assume realistic distributions of matter in the Solar neighbourhood, 
as we shall demonstrate in the next section. 

The simulations we analyse here were generated by zooming in and re- 
simulating with higher resolution a particular cluster and its surroundings formed 
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Figure 1. The left panel shows the projected mass distribution of the 
simulated halo out to the virial radius. The top panels on the right 
show the spatial distribution of particles inside a 2 kpc on a side volume 
located at 8 kpc from the halo centre, i.e. this volume is centered on 
the "Sun". The bottom panel shows their distribution in the sky. 



in a cosmological ACDM simulation, with parameters f^o = 0.3, f^A = 0.7, 
h = 0.7 and erg = 0.9. The cluster selected for resimulation is the second most 
massive cluster in the parent simulation and has a virial mass of 8.4x 10^^/i~^ Mq 
(see the left panel of Figure 1). The highest resolution resimulation of the cluster 
region has 6.6 x 10^ particles. We scale the cluster to a "Milky Way" halo by 
requiring that the maximum circular velocity is 220 km/s, which yields a scaling 
factor 7 = v^^/v^^ ~ 9.18. 



2. The phase-space structure in the Solar neighbourhood 

2.1. Spatial distribution 

One of the critical issues in understanding the outcome of the various dark- 
matter experiments consists in understanding the expected signal: Is the dis- 
tribution of particles in the vicinity of the Sun smooth or is it dominated by a 
just a few streams or even bound lumps (e.g. Moore et al. 2001)? In Figure 1 
we plot the positions of all particles inside a cubic volume of 2 kpc on a side, 
located at the position of the "Sun" . The spatial distribution of particles inside 
this representative volume is extremely smooth. This is mostly due to the fact 
that the material that ends up in the inner galaxy mostly comes from a few very 
massive halos which have rapidly mixed. 

2.2. Kinematics 

In Figure 2 we show the velocities of particles located in a box of 4 kpc on a 
side in the vicinity of the Sun. Their velocity distribution is relatively smooth, 
and appears to be quite consistent with a Gaussian (see left panel Figure 3), at 
least to the "naked eye". However if we focus on the highest energy particles 



On the phase-space structure of the Milky Way dark-matter halo 



3 



600 - 

400 - 
200 - 
- 
-200 - 
-400 - 
-600 L 




600 - 

400 - 
200 - 
- 
-200 - 
-400 - 
-600 L 






600 






400 






200 










E 



















-200 






-400 






-600 





-600-409-200 200 400 600 
V, (km/s) 



-600-409-200 200 409 600 
(km/s) 



-609499290 200 409 690 
V2 (km/s) 



Figure 2. Velocities of particles located in a box of 4 kpc on a 
side at the "Solar" radius. There are 4362 particles in this vol- 
ume, and highlighted are the 1% fastest. The velocity dispersions 
along the principal axes are (111.2, 120.1, 141.4) km/s. The lump with 
V ~ (-250, -200, -400) km/s corresponds to a halo of 1.94 x 10^° M© 
identified at z = 2.4, and accreted at z = \, or 8.2 Gyr ago. 



this seems no longer to be the shown by the particles highlighted in 

grey. The 1% fastest moving particles are strongly clumped, and their velocity 
distribution is highly anisotropic. 

To quantify the substructure present in this volume we compute the corre- 
lation function ^ in velocity space, defined as ^ = — 1, where {DD) is the 
number of pairs of particles in our simulation with velocities in a given velocity 
range (or bin), and {RR) is defined analogously for the same number of random 
points, which we draw from a trivariate Gaussian distribution determined from 
the data in the principal axes velocity frame, and is the average over ten such 
realizations. The right panel in Figure 3 shows that there is a weak excess of 
particles with similar velocities (i.e. below 100 km/s) with respect to what would 
be expected for a random multivariate Gaussian distribution. However, if we 
focus on the fastest moving particles, the excess of particles is very noticeable, 
particularly at small velocity differences, indicative of the presence of streams, 
as clearly visible in Fig. 2. 

3. Conclusions 

A dark-matter halo formed in a ACDM cosmology is not a smooth entity. Not 
only do dark-matter halos contain a large number of dark satellites, they also 
have large amounts of substructure in the form of streams. By scaling a cluster 
halo down to a galaxy halo using the ratio of their maximum circular velocities 
(equivalent to scaling their mass ratio to the one-third power), we predict that 
the Galactic dark-matter halo in the "Solar neighbourhood" should be smoothly 
distributed in space. This material is clumped in at least a few hundred kinemat- 
ically cold streams, as more detailed analyses show (Helmi, White & Springel 
2001). These streams have their origin in the different halos that merged to 
form the dark halo of the galaxy. 

Our results indicate that direct detection experiments may quite safely as- 
sume a multivariate Gaussian to represent the motions of the dark-matter par- 
ticles. However, experiments which are only sensitive to the highest energy (i.e. 



4 



Amina Helmi, Volker Springel <fe Simon D.M. White 




Figure 3. For the same box as in Figure 2 we plot the differential 
(left) and the cumulative (centre) velocity distributions. The dotted 
histograms correspond to a multivariate Gaussian. On the right, we 
plot the "correlation function" ^ as the number of neighbours with 
velocity differences in a given range compared to what is expected 
for random deviates drawn from a multivariate Gaussian. Asterisks 
correspond to ^ over all the particles in the box, whereas diamonds to 
the 1% fastest moving particles. 



fastest moving) dark- matter particles will be subject to the presence of a few 
dominating streams, and an excess of particles with similar energies should then 
be expected in comparison with a smooth Gaussian distribution. If direct de- 
tection experiments are also sensitive to the direction of motion of the incident 
dark-matter particles, the signal expected for the fastest moving particles is also 
highly anisotropic, and could be eventually be used, not just to determine the 
nature of the dark-matter, but also to recover at least partially the merging 
history of our galaxy. 
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